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O ' A partial wave analysis of the centrally produced rjir^ and r]TT~ channels has been performed 

^ . in pp collisions using an incident beam momentum of 450 GeV/c. Clear ao(980) and 02(1320) 

^ ', signals have been observed in S and D^ waves respectively. The cIPt, (p and |t| distributions 

I ' of these resonances are presented. 
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The rjTT system has been studied in several experiments in the search for exotic states [0, 0. 
The interest in this system is mainly caused by the fact that P-wave in the rjir system carries 
exotic quantum numbers J^*^ = 1 '". 

In addition to the well-known ao(980) and 02(1320), two new qq states were observed recently 
in the rjii system. In 1994 the Crystal Barrel collaboration claimed a new I'-'J^^ = 1"0^^ state, 
the ao(1450) 0, which they then confirmed in different channels [^ |^, §. This resonance could 
be considered as the isovector member of the ground state scalar multiplet instead of the ao(980), 
which is a peculiar object in meson spectroscopy and whose nature is still not completely clear, 
see for example ref. 0. In addition, the Crystal Barrel collaboration observed a new isovector 
tensor state, the 02(1660), decaying to r^vr^ g|. 

In a previous analysis of the centrally produced rjir^ system by the NA12/2 collaboration ^ 
clear ao(980) and 02(1320) signals were observed but no evidence for a J^'" = 1^^ state was 
found. An interesting feature of that analysis was the ratio of the ao(980) to a2(1320) which 
was found to be much bigger than had been observed in 7r~p interactions 0. 

This paper presents a study of the centrally produced rjir final state. Firstly the reaction 

pp ^ Ps{r]7T°)pf (1) 

has been studied at 450 GeV/c, where the r] is observed decaying to 27s . The subscripts / 
and s indicate the fastest and slowest particles in the laboratory respectively. The WA102 
experiment has been performed using the CERN Omega Spectrometer, the layout of which is 
described in ref. [1^. The GAMS-4000 multiphoton spectrometer has been used to detect 7s . 



Reaction (|T]) has been isolated from the sample of events having two outgoing charged tracks 
and four 7s by first imposing the following cuts on the components of the missing momentum: 
|missing Pj.\ < 17.0 GeV/c, |missing Py\ < 0.16 GeV/c and |missing P^l < 0.12 GeV/c, where 
the X axis is along the beam direction. A correlation between pulse-height and momentum 
obtained from a system of scintillation counters was used to ensure that the slow particle was 
a proton. 

The separation of events belonging to reaction (|l|) has been performed on the basis of a 
kinematical analysis (6C fit, four-momentum conservation being used and the masses of the two 
mesons being fixed). There is evidence for a small A"*" (1232) signal in the p/vr" mass spectrum 
which has been removed by requiring M (p/vr^) > 1.3 GeV. Detector acceptances and efficiencies 
have been calculated using a Monte-Carlo simulation taking into account the geometry of the 
set up, detector resolutions, event selection criteria and the kinematical fitting procedure. The 
resulting centrally produced rjir^ effective mass distribution corrected for efficiency and rescaled 
to the total number of observed events is shown in fig. |I]a) and consists of 6045 events. There 
is clear evidence for the ao(980) and 02(1320). As can be seen from the 77 mass spectrum 
shown as an inset (fig |l]b)) there is approximately 30 % background below the 7] signal. To 
investigate the effects of this background the side bands around the t] signal have been studied. 
Superimposed on the mass spectrum as a shaded histogram is the estimate of the background. 
Fig. ^) shows the background subtracted mass spectrum. 

To determine the parameters of the resonances produced in reaction (|1|) a fit to the efficiency 
corrected mass spectrum has been performed using a parametrisation of the form 

dN/dm (m) = BG{m) + ai\BWa,i9so)im,0)\'^ + a2\BWa,{i320)im, 2)1'^ 

2 



where 

BG{m) = (m - rrithrTe-'^"'-^"'' 

represents the background, m is the rjii mass, nithr is the rjn threshold mass and a„, a, (3 and 7 
are parameters to be determined from the fit. Relativistic Breit-Wigner functions BWi{m, J), 
where J is the spin of the resonance, have been convoluted with a Gaussian to account for 
the experimental mass resolution. The Breit-Wigner used to describe the ao(980) uses the 
Flatte formula [^. In order to describe the centrally produced rjir mass spectra the differential 
cross section has been multiphed by the kinematical factor (M^jn — Amlf^^Y^'^ /M^^ |1^. The 
fit is shown in fig. |l]c) and gives M(ao(980)) = 975 ± 7 MeV, r(ao(980)) = 72 ± 16 MeV and 
M(a2(1320) = 1308 ± 9 MeV, r(a2(1320)) = 115 ± 20 MeV. The parameters found for the 
ao(980) and 02(1320) are consistent with those from the PDG [0. 

A PWA has been performed in the mass interval from 670 to 2000 MeV for the rjir^ final state 
in 60 MeV mass bins. The PWA of the centrally produced rjir system has been made assuming 
that the t^tt system is produced by the collision of two particles (referred to as exchanged 
particles) emitted by the scattered protons. The z axis is defined by the momentum vector of 
the exchanged particle with the greatest four-momentum transferred in the rjir centre of mass. 
The y axis is defined by the cross product of the momentum vectors of the two exchanged 
particles in the pp centre of mass. The two variables needed to specify the decay process were 
taken as the polar and azimuthal angles (^,0)=f2 of the tt in the r/vr centre of mass relative 
to the coordinate system described above. The amplitudes used for the PWA are defined in 
the refiectivity basis [Q. In this basis the angular distribution /(fi) is given by a sum of two 
non-interfering terms corresponding to negative and positive values of refiectivity e. The waves 
used were of the form Jq (me = 0—), J_ (me = 1—) and J+ (me = 1-I-). Only S", P and D 
waves (corresponding to J = 0, 1,2) have been used in the PWA with ?7i = 0, 1. It has been 
determined from a study of the moments that the contributions of higher waves are negligible 
in the mass region under study. Finally, 12 parameters have to be determined from the fit 
to the angular distributions. In this model of the PWA the number of non-trivial solutions 
in each mass interval can be less than or equal to 8 due to the ambiguities problem W^ 11^ 



In each mass bin an event-by-event maximum Likelihood method has been used. The PWA 
analysis has been performed by extending the likelihood function to include the background 
subtraction, namely the function 

F = -J2 ln{Iin)} + J2 tLMeiM + E in{i{n)} 

i=l L,M i=l 

has been minimised, where A^ is the number of events in a given mass bin, eiM are the efficiency 
corrections calculated in the centre of the bin, tLM are the moments of the angular distribution 
and Nfjg is the number of background events. 

All solutions for S, P and D waves have been found for each interval. In the mass region 
above 1.2 GeV, because the D-wave is the dominant contribution, the solutions are not affected 
by ambiguities and the 8 solutions are effectively identical. The threshold region does suffer 
from ambiguities, however, we have picked the solution in which the known ao(980) is in the 
S-wave. 

The physical solution for the rj-K^ system is shown in fig. [T|d). As can be seen the D^ wave 
dominates above 1.2 GeV. The fit of the D^ wave amplitude squared with a Breit-Wigner 



function gives parameters for 02(1320) meson similar to the ones from the fit to the efficiency 
corrected mass spectrum. A fit to the S wave amphtude squared gives parameters for ao(980) 
similar to those from a fit to the mass spectrum. There is no evidence for an ao(1450) nor an 
02(1660). As can be seen there is no evidence for other significant structures in any other wave. 

In previous analyses we have observed that 1 = 1 mesons are suppressed in the reaction 
pp -^ ppX^ but are enhanced in the reaction pp — > A'^^pX" [^. Therefore, the reaction 



pp -y ps{r]7r )A++ (2) 

with the A^"*" observed decaying to PfT^~^, has been studied at 450 GeV/c. The t] is observed 
decaying to 27s . Reaction (g) has been isolated from the sample of events having four outgoing 
charged tracks and two 7s by imposing the cuts on the components of missing momentum used 
for reaction (|l]). Events have been separated by a kinematical analysis (5C fit). The events 
with a A^+(1232) in the final state have been selected by requiring M (pjvr^) < 1.4 GeV. The 
background process pp —>■ Ps{'ri'K^'K^)pf gives the biggest background below the A++ peak. By 
requiring M (r^Tr+vr^) > 1.5 GeV the signal to background in the A++ region is greater than 10 
(figlpD)). The resulting acceptance corrected 777?" mass spectrum rescaled to the total number 
of observed events is shown in fig. ^p) and consists of 8027 events. The mass spectrum is similar 
to that observed in the t^tt" final state but the ratio ao/a2 is different. In the 77 mass spectrum 
(not shown) there is approximately 30 % background below the r] signal. To investigate the 
effects of this background the side bands around the rj signal have been studied. Superimposed 
on the mass spectrum as a shaded histogram is the estimate of the background. 

A fit to the centrally produced background subtracted 777?" system is shown in fig|^c). The 
fit yields M(ao(980)) = 988 ± 8 MeV, r(ao(980)) = 61 ± 19 MeV and M(a2(1320) = 1316 ± 9 
MeV, 1(02(1320)) = 112 ± 14 MeV. 

A PWA has been performed in the mass interval from 670 to 2000 MeV for the r/vr" final 
state in 60 MeV mass bins using the method described above. The physical solution for the 
?77r~ system is shown in fig. ^d). As can be seen the D+ wave dominates above 1.2 GeV, the 
fit of the D+ wave amplitude squared with a Breit-Wigner function gives the parameters for 
02(1320) meson similar to the ones from the fit to the efficiency corrected mass spectrum. A fit 
to the S wave amplitude squared gives the same parameters for ao(980) similar to those from 
the fit to the mass spectrum. Again there is no evidence for an ao(1450) nor an a2(1660). As 
can be seen there is no evidence for other significant structures in any other wave. 

The fits to rj-K^ and r/Tr" mass spectra, give the following values for the production ratio of 
the ao(980) and a2(1320) decaying to tjti : 

a{pp -^ pp[ao(980) -^ tj-k]) 



for reaction (|l]) and 



a{j)p -^ pp[a2(1320) — > rj-n]) 
a{pp -^ pA++[ao(980) -^ r/vr]) 



2.0 ±0.3 



0.8 ±0.2 



a{pp -^ pA++[a2(1320) -^ r/vr]) 

for reaction (^. In the charge exchange reaction (7r~p — > rjTf^n), where the 02 (1320) meson 
dominates the rjir^ mass spectrum [Q: 

a{-n~p -^ [ao(980) -^ rjTx]n) 



cr(7r p -^ [02(1320) -^ rjTi]n) 
4 



0.15 



at 38 GeV/c beam inomentum. 

In previous WA102 analyses it has been observed that centrally produced states have dif- 
ferent (IPt dependencies, where dPT is the difference in the transverse momentum vectors of 
the two exchange particles [|l^]. The ratio of the production cross section for dPx < 0.2 GeV 
to cIPt > 0.5 GeV is significantly different for qq states and for the glueball candidates. It has 
been observed that all undisputed qq states which can be produced by DPE have very small 
values for this ratio (< 0.1). States which cannot be produced by DPE (with negative G parity, 
for example) have slightly higher values (~ 0.25). And all non-gg candidates /o(980), /o(1500), 
/o(1710) and /2(1950) have values of this ratio about 1. A study of the r/vr^ and r/vr" systems, 
which can not be produced by DPE, has been made as a function of (IPt- The results are given 
in table 0. 

In addition, an interesting effect has been observed in the azimuthal angle which is defined 



as the angle between the px vectors of the two outgoing protons [|16|. In order to determine 
the azimuthal angle for the ao(980) and a2(1320) the tjtt^ and r/vr" mass spectra have been 
fitted in 30 degree bins of (p with the parameters of the resonances fixed to those obtained from 
the fit to the total data. The resulting distributions for the aQ(980) and a2(1320) are shown in 
fig. |]a) and b) respectively. The distributions for the ao"(980) and 02^(1320) are shown in fig. || 
a) and b) respectively. 

In order to determine the four momentum transfer dependence (t) of the resonances ob- 
served, the ?77r° mass spectrum has been fitted in 0.1 GeV^ bins of t with the parameters of the 
resonances fixed to those obtained from the fits to the total data. Fig. ^) and d) shows the four 
momentum transfer from one of the proton vertices for the ao(980) and 02(1320) respectively. 
The distributions have been fitted with a single exponential of the form exp{—b\t\) and the 
values of b found are 6 = 6.2 ± 0.8 for the aj](980) and 6 = 8.8 ± 0.4 for the a^(1320). 

In order to determine the four momentum transfer (|t|) at the beam-A"*""^ vertex for the 
0(7(980) and a^(1320) the rjir' mass spectrum has been fitted as described above. The resulting 
distributions are shown in fig. |c) and d) for the a^(980) and a^(1320) respectively. The 
distributions have been fitted with a single exponential of the form exp{—b\t\) and the values 
of b found are 6 = 4.0 ± 0.3 for the 0(^(980) and 6 = 5.9 ± 0.5 for the a^(1320). 

The four momentum transfer (|ts/oto|) distributions at the target-slow vertex are shown in 
fig. |e) and f) for the ao{980) and a^(1320) respectively. The data for |t| < 0.1 GeV^ has 
been excluded from the fit due to the poor acceptance for the slow proton in this range. The 
distributions have been fitted with a single exponential of the form exp{—b\t\) and the values 
of b found are 6 = 7.5 ± 0.4 for the 0(^(980) and b = 7.9± 0.6 for the a^(1320). 

In summary, a partial wave analysis of the centrally produced rjir^ and rjir^ channels has 
been performed in the mass region 670 — 2000 MeV. Clear ao(980) and 02(1320) signals are 
seen in the S and D+ waves respectively. No evidence is found for any significant structures in 
other waves. In particular, there is no evidence for the ao(1450) nor the a2(1660). The ratio of 
the ao(980) to 02 (1320) is found to depend on the production mechanism. 
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Table 1: Production of the resonances as a function of dP'^ expressed as a percentage of their 
total contribution and the ratio (R) of events produced at dPr < 0.2 GeV to the events produced 
at dPT > 0.5 GeV. 



Resonance 


dPT<0.2 GeV 


0.2<dPT<0.5 GeV 


dPT>0.5 GeV 


p dPT<0.2GeV 
-"- dPT>0.5GeV 


a[](980) 


25±3 


33±5 


42±4 


0.57 ±0.09 


a^(980) 


14±3 


37±2 


49±2 


0.29 ±0.06 


a^(1320) 


10±2 


38±2 


52±3 


0.19 ±0.04 


a^(1320) 


9±3 


39±2 


52±2 


0.17 ±0.06 



Figures 

Figure 1: a) The efficiency corrected mass spectrum for the centrally produced rj-K^ events. 
Superimposed as a shaded histogram is an estimation of the background contribution. Inset b) 
the 77 mass spectrum, c) The efficiency corrected background subtracted mass spectrum for 
the centrally produced ?77r° events. The curve is the result of the fit described in the text, d) 
The physical solutions from the PWA of the rjTC^ events. 



Figure 2: a) The efficiency corrected mass spectra of centrally produced rjir^ events. Inset b) 
presents the distribution to the p/7r+ invariant masses, c) The efficiency corrected background 
subtracted mass spectrum for the centrally produced rj-K" events. The curve is the result of the 
fit described in the text, d) The physical solutions from the PWA of the r^vr" events. 



Figure 3: The azimuthal angle between the pt vectors of the slow and fast protons for a) 
the aQ(980) and b) the 02(1320). The four momentum transfer distributions for c) the aQ(980) 
and d) the 02(1320), with the fits to the form e"**'*'. 



Figure 4: The azimuthal angle between the pt vectors of the slow proton and fast A"^"*" 
for a) the a^(980) and b) the a^(1320). The four momentum transfer distributions from the 
beam-fast vertex for c) the aQ"(980) and d) the 02^(1320) and for the target-slow vertex e) the 
0(^(980) and f) the a^(1320), with the fits to the form e~^l*l. 
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